THEORY OF SUSPENSION SEGREGATION IN PARTIALLY 
FILLED HORIZONTAL ROTATING CYLINDERS 

Rama Govmdarajan, 1 '^ Prabhu R. Nott, 2 'Q and Sriram Ramaswamy* 3 '! 

11 Fluid Dynamics Unit, Jawaharlal Nehru Centre for Advanced 
Scientific Research, Jakkur, Bangalore 560 064, India 

2 Department of Chemical Engineering, 
Indian Institute of Science, Bangalore 560 012, India 

3 Centre for Condensed- Matter Theory, Department of Physics, 
Indian Institute of Science, Bangalore 560 012, India 

Abstract 

It is shown that a suspension of particles in a partially-filled, horizontal, rotating cylinder is lin- 
early unstable towards axial segregation and an undulation of the free surface at large enough par- 
ticle concentrations. Relying on the shear-induced diffusion of particles, concentration-dependent 
viscosity, and the existence of a free surface, our theory provides an explanation of the experiments 
of Tirumkudulu et al, Phys. Fluids 11, 507-509 (1999); ibid. 12, 1615 (2000). 
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The primary motivation for this work is to understand the experiments of JI], 0], showing 
spontaneous segregation in sheared suspensions. In these experiments, monodisperse poly- 
mer spheres several hundred microns in diameter were suspended uniformly in fluids about 
1000 times as viscous as water, with the same mass density as the particles, so that there 
was no sedimentation. The initial experiments [[I]] were carried out with a suspension in a 
horizontal Couette cell, i.e., in the gap between two concentric cylinders, but in subsequent 
studies H, a single horizontal cylinder of radius R was used. In both sets of experiments 
the container was filled only partially, i.e., there was a free surface. Let us restrict our 
description to the experiments in for simplicity. When the cylinder was rotated at a 
tangential speed Vq = RQ about its symmetry axis the initially uniform suspension was 
found to undergo a dramatic instability towards segregation into bands of high and low 
concentration, with wavevector along the cylinder axis. The surface profile was modulated 
as well, i.e., the thickness of the fluid layer varied along the axis (see Fig. p]), with thicker 
regions corresponding to higher concentration. No instability was seen if there was no free 
surface, i.e. when the cylinder was completely filled with suspension. 

The phenomena reported in the experiments Jl|, Q are among the many intriguing effects 
known to arise in suspensions of non-Brownian particles in highly viscous fluids, driven 
by shear flow or sedimentation. The feature of these suspensions that is of relevance in 
this paper is that the particles diffuse |J even though their thermal Brownian motion is 
negligible. The microscopic explanation |3], ^, |5| for this diffusion is that the hydrodynamic 
interaction between the particles renders their motion chaotic, even in the Stokesian limit 
where inertia of fluid and particles is ignored. The diffusive flux of particles has two parts, 
one driven by a gradient in the particle volume fraction 0, the other by a gradient in the 
deformation-rate 7. For situations in which the predominant variation is with respect to 
a single coordinate z (e.g., the axial coordinate of the cylinder in [l], ||) and time t the 
local volume fraction <j)(z, t) (integrated over the remaining directions) of particles obeys the 
conservation law 

dt dz' 1 ' 

where the shear-induced current || || can be written as 

j = -/ c (0)a 2 7 ^-/ s (0)a 2 0|j, (2) 

In (Q), a is the particle radius, and f c and f s are dimensionless functions of the particle 
volume fraction [H, EL f7|. Note that (0) says that particles can move in the absence 
of concentration gradients, or even against concentration gradients, if the gradient in the 
deformation-rate is appropriately directed. These equations are an essential ingredient of 
our theory of the shear- induced segregation seen in [fll 0. 

The main result of our analysis is that equations ([!]) and (@), when applied to neutrally 
buoyant Stokesian suspensions in horizontal rotating cylinders, predict precisely the insta- 
bility seen in the experiments of if the concentration is large enough. The growth rate 
T q of the instability varies as q 2 for small wavenumber q. The parameters which govern the 
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instability depend only on the volume fraction of the suspension and the fill fraction of 
the cylinder. At low rotation rates, Q scales completely out of the problem: the range of 
unstable wavenumbers is therefore independent of Q. For just above the instability-onset 



value C , T q reaches a maximum at q = q* ~ y (0 — 4>c)pg/o~ where cr is the surface tension 
of the suspension. For reasonable values of these parameters we find that the fastest growing 
mode has wavelength of order centimetres. We also explain why the instability disappears 
for large rotation rates. 

We now obtain coupled equations of motion for the particle concentration and free-surface 
profile, and show that these lead naturally to the above results. Consider a homogeneous 
suspension with kinematic viscosity u, filling a fraction F of the volume of a horizontal 
cylinder of radius R, rotating about its symmetry axis with a tangential velocity Vq = 
RQ (See Fig. 0). Recall first the results of ||. (a) The dimensionless combination (3 = 
F^JgR 2 / uvq, where g is the acceleration due to gravity, measures the relative importance 
of gravitational and viscous forces; (b) As vq is increased (i.e., as (3 is decreased), a fluid 
film of thickness w is dragged up and coats the cylinder wall; (c) Since w, for low speeds, 
is smaller than the depth of the residual pool of fluid at the bottom of the cylinder, the 
thickness profile has a "bump" at the bottom, (d) Once w reaches a value ~ FR, which 
occurs for (3 = (3 C ~ 1.4, the growth of w saturates since all the available fluid then coats 
the cylinder more or less uniformly, and the bump disappears, (e) For higher speeds, i.e., 
for (3 < (3 C , w is effectively independent of the rotation speed, and is determined simply by 
the geometrical statement w/R = F. The mechanism we propose below for the instability 
applies only when the thickness is determined by the rotation speed by an explicit balance 
between viscous and gravitational forces which is why, in 0, the instability disappears when 
the bump does. 

Consider a general situation (Fig. [|) where the thickness w(z, t) of the fluid film dragged 
up, as well as the volume fraction field <p(z, t) (and hence the viscosity), are varying in space 
and time. The component of the deformation rate that could vary in the axial direction is 
given by the velocity difference across the layer divided by the thickness: 



where «2 is a pure number of order unity, independent of material parameters. The experi- 
ments of JL|, 0] are performed on highly viscous fluids, so that the Reynolds number is very 
small over the entire range of speeds and length scales studied. We shall therefore work 
in the limit of zero Reynolds number, where the inertia of particles and fluid are ignored. 
Accordingly, the balance of gravitational, viscous and interfacial forces per unit area of the 
layer tells us that 





(4) 



and in particular that the layer thickness in the steady, spatially uniform state is 




(5) 
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and (|5|), p, r]((j)), and a are respectively the density, effective viscosity (as a function 
TOfl of the local particle volume fraction 0) and surface tension of the suspension, g is 



the acceleration due to gravity, and ai is another geometrical factor of order unity, with 
no dependence on any material parameter. The fill fraction F determines the angle made 
by the cylinder with the free surface of the pool of suspension, and hence the details of the 
flow in the fluid layer. This is reflected in the parameters a.\ and a.2 in our model, but our 
conclusions are qualitatively insensitive to their precise numerical values. 

Let us now perturb the thickness and concentration fields: [w(z,t), 4>(z,t)] = [w + 
Sw(z, t), 0o + VK- 2 ; t)}- This will in turn lead to perturbations of the local values of 7 and rj, 
yielding closed equations of motion for the evolution of Sw(z,t) and ip(z,t) via flU) and (§]). 
We work to linear order in 5w and ip. Let us work in terms of the nondimensional quantities 
H = 5w/w, t = (v /a2w)t, ( = z/a. -0 is of course already dimensionless. Note that to 
write (HD and (§) in terms of thickness and concentration fluctuations, we must use @ to 
express the local deformation-rate in (fj) in terms of the thickness perturbation, and replace 
a local viscosity perturbation by a local concentration fluctuation via 6r)/r) ~ J\fip where 

d In 77 



(6) 



Although the procedure is straightforward, some care must be taken in obtaining the per- 
turbation equation from (^J): the perturbation of 7 will involve w, which must then once 
again be eliminated in favour of 5w, ip. Carrying out these steps, and Fourier-transforming 
with respect to (, i.e., considering spatial variation of the form expiq(, we find that the 
Fourier components H q , %fj q obey 
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«20o/s(0o), C = a 2 / c (0o) and S = -. (8) 

pga 2 

The stability or otherwise of our sheared suspension is determined by the characteristic 
equation 

A 2 + (2 + Dq 2 )\ + Eq 2 + CSg 4 = (9) 
for the eigenvalues A of the dynamical matrix M in (0), where 

D = £ + C - US, 

2C - ATS. (10) 



E 



For q — > 0, the solutions to (P) are 



Ai 



E 2 

y 2 



D)q 2 . 



(11) 
(12) 
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We see from (|Tl| ) that the uniform state is linearly unstable to segregation and thickness 
modulation if E < 0, with perturbations growing at a rate T q ~ q 2 at small q, and particles 
concentrating in the thick regions. It is straightforward to show that Ai turns over at larger 



q, passing through zero at q = J—E/CY>, with a peak located, for E — > , at 



q* ~ J-E/2CE (13) 



> 2 J^tA (14) 



which determines the observed wavelength of the initial instability |IlJ. For a given fill 
fraction, E can vary only with the volume fraction 0. If there is an instability, it must then 
be because E turns negative, in general as C — 0, as crosses a critical value <f) c . This leads 
to the main result presented at the start of the paper. In terms of the parameters in (0) and 
(§), the instability criterion is 

dr] 
i] d<p 

which should in general happen in real suspensions at large enough O . 

In more detail, note that the coefficient / c (0) originates || [| from a direct shear-induced 
self-diffusion as well as from a tendency to move down viscosity gradients. The latter ten- 
dency opposes the instability, as we shall now show. For a Newtonian suspension, rj varies 
only if does. Thus, by the arguments of Leighton and Acrivos 0], we can write the current 
in (|I|) as 

j = -a 2 0[M s 0^ + 7(M C + M v H)^i (15) 

where M c , M s and M v are order-unity phenomenological quantities. Comparing ( |T5| ) and 
(0) we see that f s = 4>M S and f c = (p(M c + J\fM v ), and the instability criterion (|TJ|) becomes 
(M s — 2M rj )M > 2M C . The experiments of H § have determined these coefficients for 
diffusion in the gradient direction, while diffusion in the problem we consider is in the 
vorticity direction. Further, a microscopic theory for determining them is also not available. 
Provided M s > 2M V , the growth of the viscosity with |], 1(| means that the instability 



should always arise at large enough concentration. An independent measurement of these 
coefficients is clearly called for. 

We now assume that the uniform state is unstable and ask whether the typical wavenum- 
ber of the segregation changes from its initial value g*, and whether nonlinear terms cause 
the exponential growth to saturate at long times. It suffices to look at the dynamics for the 
slow mode (eigenvalue Ai in the linearised limit). In that case H q in (0) is "slaved" to ip q , 
H q ~ —Afip q /(2 + Sg 2 ). Expanding (|2]) about 0o ; an d retaining terms nonlinear in ip q , we 
obtain the effective equation of motion 

di) d 2 ,E , , 2 (3 , NS&rp 



where C2 and C3 are coefficients arising from the dependence of f c , f s , and rj. Eqn. ( |16|) 
is well-known in the domain-growth literature (12], [13] . In particular, it has been shown |14[] 
that the nonlinear terms in flTEp with E < cause the characteristic wavelength L(t) of the 
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pattern of segregation at time t to grow as hit at long times. This extremely slow growth 
should in principle be testable by a patient experimenter. 

We have shown that two classic properties of non-Brownian suspensions, viz., 
concentration-dependent viscosity || [HJ and shear-induced diffusion [[|, lead to a natu- 
ral explanation of the experiments of [|TJ on segregation in suspensions in partially filled, 
rotating horizontal cylinders. Our dynamical equations, at high enough concentration, dis- 
play an instability towards axial segregation and a modulation of the free surface, with 
particles accumulating under the crests of the modulation. For parameter values, say, 10 % 
past the instability, taking a plausible surface tension of 30 dyne/cm, implies a wave- 
length of about 3 cm for the fastest growing mode at onset, which is consistent with the 



experiments of 0. Independent measurements of the parameters in (|Tj) and (0), in tran- 
sient experiments, should provide a stringent test of our theory, as should studies of the 
long-time behaviour of the wavelength of the segregation pattern. 
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FIGURE CAPTIONS 

Figure 1: Schematic of steady-state surface profile of a horizontal cylinder of radius R, 
rotating with angular speed Q. The thickness w of the suspension layer dragged up, as well 
the concentration of solute indicated by shading, are modulated with respect to the axial 
coordinate z (see ref. 0). 

Figure 2: Cross-sectional schematic of profile of fluid layer of characteristic thickness w 
and kinematic viscosity u, dragged up against gravity g by the wall of a cylinder, rising at 
tangential speed Vq. 
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FIG. 2: 
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